A new set of variables in the three-body problem 
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We propose a set of variables of the general three-body problem both for two- 
dimensional and three-dimensional cases. Variables are (A, 9, A, O, k, to) or equiva- 
lcntly (A, 9, L, I, k, o>) for the two-dimensional problem, and (A, 9, L, I , k, u>, 4>, ip) for 
the three-dimensional problem. Here (A, 6) and (A, 8) specifies the positions in the 
shape spheres in the configuration and momentum spaces, k is the virial ratio, L is 
the total angular momentum, / is the time derivative of the moment of inertia, and 
u>, <p, and ip arc the Eulcr angles to bring the momentum triangle from the nominal 
position to a given position. This set of variables defines a shape space of the three- 
body problem. This is also used as an initial condition space. The initial condition 
of the so-called free-fall three-body problem is (A, 6, k = 0, L = 0, 1 = 0, w = 0). We 
show that the hyper-surface / = is a global surface of section. 
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1 Introduction and motiva- 
tion 

The three-body problem has been studied exten- 
sively by numerical simulations and by analytical 
methods. However, in full generality this problem 
is of too high dimensionality and is too complicated 
for a systematic analysis. This has led so far to 
the study of various simplified and restricted ver- 
sions of the problem. The free-fall problem (|T], 
[T4] . referred to as the FFP; the definition will be 
given later) is one of these. The isosceles prob- 
lem ([2], [S], [H]) and the rectilinear problem ([5], 
[151 f!6] ) are other examples. A research looking 
at the whole phase space of the planar or three- 
dimensional three-body problem is desirable. 

In the present work, we propose a new setting 
of the problem which is suitable for large-scale nu- 
merical studies, and which is hopefully suitable for 
theoretical consideration. Our setting is summa- 
rized into a word 'shape space'. The shape space is 
a direct product of the shape sphere in the config- 
uration space, the shape sphere in the momentum 
space, and their relative size and orientation. This 



is not the phase space, but can be expanded to re- 
cover the phase space. 

In getting the shape space, we have two guiding 
principles: equivalence relation and boundedness. 
We do not want to integrate orbits which can be 
transformed into one another by a suitable change 
of variables. These orbits are considered equivalent. 
We will consider only orbits of different equivalence 
classes. Orbits belonging to different equivalence 
classes will be called independent. Boundedness is 
related to the size of the initial condition space. In 
order that a numerical study of the totality of the 
solution of the three-body problem be feasible, the 
initial condition space must be either bounded or 
finite. Otherwise, it is impossible to exhaust the 
initial conditions using even the fastest computer. 
So we would like to impose boundedness on the 
ranges of variables. 

We consider the planar three-body problem since 
the extension to the three dimensions is not so dif- 
ficult. The phase space is R 12 . As is well-known 
(|18j, p. 351), the possible minimum order of differ- 
ential equations is four. In order to attain this or- 
der, several processes arc consecutively done. Thus, 
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equations of motion are reduced from the 12th to 
the 8th order by using the four integrals of motion 
of the center of gravity. The use of the angular 
momentums reduces the order to 7, and use of the 
elimination of the nodes reduces the order to 6. 
Lastly, it is possible again to reduce the order of 
the equations to four by using the integral of en- 
ergy and eliminating the time. 

This conventional reduction, however, does not 
lead to a best choice of variables for the initial 
value problem. The space of the states of motion 
described by these 4th-ordcr differential equations 
may have extremely complicated structure. Even 
if we stop the reduction at the 6th-ordcr, variables 
have infinite ranges. 

The concept of the shape sphere in the configura- 
tion space seems due to McGehee in the rectilinear 
three-body problem when he devised the variables 
now called after his name [H]. Later Moeckel [I3j 
explicitly introdcued the shape sphere in the pla- 
nar three-body problem. Then, this notion is essen- 
tially used in [7] in order to obtain the figure-eight 
periodic solution with variational techniques. Our 
shape space extends the notion to the whole phase 
space. 

We start in §2 by introducing the FFP and the 
involved study for it. In §3 we extend the FFP to 
any given mass ratio and the whole phase space in 
the planar case. In §4 we suggest the existence of 
the 'semi' global surface of section and generalize 
this results to the three dimensions. In the final 
section, we summarize and discuss our results. 

2 The Free-Fail Problem 
2.1 Definition 

The free-fall problem (FFP) is characterized by 
the zero initial velocities, and has been extensively 
studied by Russian and Japanese schools 
[Hj. [17]), In the FFP, the total energy of the 
three bodies m,i,i = 1,2,3 is negative and their 
angular momentum is zero. We here consider the 
equal mass case: mi = 777,2 = 777.3 ■ I n this problem, 
motions starting from similar triangles transform 
into one another under appropriate changes of co- 
ordinates and time, so we identify these motions. 
Dissimilar triangles correspond to independent mo- 
tions. 



Let mass points 7772 and 777,3 stand still at 
A(-0.5,0.0) and £(+0.5,0.0), respectively in the 
(x, y) plane and mi stand still at a point P{x, y) 
where 

{x,y)eV= (1) 
{{x, y) : x > 0, y > 0, (x + 0.5) 2 + y 2 < 1}. 

If mi changes position in T>, then triangles satis- 
fying the condition 7772777,3 > m\mi > 7771777,3 are 
exhausted. Conversely, any triangle is similar to 
one of the triangles formed by three mass points 
mi, 777,2, and 777,3 as above. Thus the positions of 
P G T> specify all possible initial conditions. 

2.2 An attempt to include velocities 

Anosova et al. ([3]) considered the region T> of the 
free-fall problem. Then they supposed that the sys- 
tem rotates in the plane of initial triangle (2D prob- 
lem) counterclockwise; the velocity vectors of com- 
ponents A (distant component) and C (component 
inside T>) arc orthogonal to their radius-vectors in 
the center-of-gravity coordinate system; the angu- 
lar momenta of these bodies are the same; the ve- 
locity of component B is given so that the center-of- 
gravity of the triple system is motionless; the speed 
of rotation is parameterized by initial virial ratio k. 
Thus the initial conditions are defined by three pa- 
rameters: coordinates (x, y) of C component in the 
region T> and virial ratio k. 

However, their formulation lost the boundedncss 
of the initial configuration space. This boundedncss 
is one of the most important properties of the FFP. 
What wc should do is to recover this. 

3 The definition of our vari- 
ables for the planer case 

3.1 Equations of motions for the pla- 
nar three-body problem 

Let mfc > be the masses of point particles 
with positions G R 2 and momenta pk G 
R 2 ;fc = 1,2,3. Let q, p G R 6 denote the vectors 
(qi,q2,q3), (Pi,P2,P3)- The three-body problem 
is governed by the Hamiltonian function 

H(p, q) = |p ■ A^p - U(q) = T(p) - U(q) (2) 
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where A is the 6x6 mass matrix 
dia<7(mi, mi, 77Z2, TO2, m3, TO3), a dot denotes 
the scalar product in R 6 , and 

TO1TO2 m 2 m 3 m z mi 

U{q) = -, |+i r + i 1- (3) 

|qi-Q2| Iqa-qal |<I3 — qi| 

Hamilton's equations are 

q = A" 1 p, p = Vtf(p). (4) 

where a dot above letters denote the time deriva- 
tive. Without loss of generality, we can assume the 
center of gravity remains at the origin: 

3 3 

rriiCii = 0, ^ pi = 0. (5) 

i=l i=l 

All these four equations characterize the planar 
three-body problem. 

3.2 Reduction to our variables 

The dimension of the original phase space is twelve. 
The restriction ([5|) reduces it to eight. The four 
variables out of eight are for the configuration 
space, and the remaining four are for the momen- 
tum space. The restriction is equivalent with 
the fact that the two sets of the three vectors form 
two triangles. In the configuration space, there re- 
main two variables. Obviously, these are the size 
and orientation of a triangle. Here, the orientation 
means the direction angle of a selected edge with 
respect to the coordinate axis. In the momentum 
space, there also remain two variables. These are 
the size and orientation of a triangle. 

Let us now look for the dimension of the space 
in which all the independent orbits are contained. 
Here, we say that two orbits are dependent (resp. 
independent) if they can (resp. can not) transform 
to each other under coordinate and time transfor- 
mations. Let us express the four variables in the 
configuration by (f, r, to), where and f is a two di- 
mensional vector and represents the form of the 
triangle, r is the size, and to is the orientation. In 
a similar manner, we express the four variables in 
the momentum space by (F,i?, f2). Thus we have 
eight variables (f , F, r, R, to, Q). 

Now, consider two sets of variables 
(f , F, r, R, a, O) and (f, F, r', R', to, Ct). If r' = ar 



and R' = (3R for appropriate constants a and (3 
(this represents a scale transformation), motions 
starting at the two initial conditions are not 
independent. If we choose a particular transfor- 
mation r, R — ► then the variables reduce to 
(t,F,l,R*,to,Q), that is, (f,F,R*,to,Cl). 

Two orientations to and are not independent. 
In fact, to is measured with respect to a fixed direc- 
tion in the configuration space, and Q is measured 
with respect to a fixed direction in the momen- 
tum space. However, the axes in the momentum 
space can be adjusted to those in the configura- 
tion space. Let us consider two sets of variables 
(f,F,R*,to,(l) and (f,F,R*,to',Cl'). Then, as is 
easily understood, if f2 — to = f2' — to' (rotation), 
motions starting at the two initial conditions are 
not independent. If we choose a particular rotation 
to, f2 — > 0, r2* , variables reduce to (f , F, R*, 0, fi*), 
that is, (f , F, R*, Q*). This is our variables. 

In the following two subsections, we will carry 
out the above program to the final set of variables. 

3.3 From Shape Plane to Shape 
Sphere 

The variables in the position space in this subsec- 
tion is equivalent with those for the shape sphere 
[13] [7] and for the FFP [I]. In [Hj, the shape 
sphere has been defined as the sphere of constant 
moment of inertia originally introduced by McGe- 
hee [5]. Our definition is slightly different from it. 
Our construction of the shape sphere is naive and 
intuitive. 

We here connect the representations of triangles 
on the plane and on the sphere. In order for this, 
let us introduce the shape plane as follows. Let 
mass points 7712 and m.3 be at B(— 0.5,0.0) and 
C(0.5,0.0), respectively in the (x,y) plane and the 
position of mi be P(x, y). If mi changes position in 
the plane, then triangles exhaust all of the shape. 
We call this the shape plane. 

Obviously the points A(0,&) and A'(0,-^) 
correspond to the equilateral configurations. The 
x-axis corresponds to collinear configurations, 
whereas the y-axis to the isosceles triangles in 
which the length of two edges m\mi and mim^ 
are the same. If mi is on B, C or its distance from 
the remaining two tends infinity, then the shape 
corresponds to binary collision £3, Ei or E\ where 



3 



Ei means that particle i goes to infinity or parti- 
cles (j, k) collide where i ^ j and i ^ k. If three 
masses are equal and mi is at D{— 3/2,0), 0(0,0) 
or 25(3/2, 0), then the corresponding configurations 
is C2 , C\ or C3 where Ci represents the isosceles or 
collincar configuration. 




Figure 1: The relation between shape plane and 
shape sphere. 

Now, let us obtain the transformation between 
(x,y) on the shape plane and (A, 8), the longitude 
and latitude, on the shape sphere. We put the 
sphere of radius in the three-dimensional space 
(x, y, z) with center at N(0, 0, as in Fig. [TJ The 
equation for the shape sphere is 



y 2 (z-Vz/Af 

(V3/4) 2 (V3/4) a (V3/4) 2 



1. 



Every straight line connecting 2?i(0, 0, and 
(cc,y,0) meets the sphere at a point. Hence every 
point on the shape plane is mapped to a point on 
the shape sphere, and vice versa. (As usual, infinity 
in the (x, y) plane is treated as a point.) We denote 
several particular points on the shape sphere as fol- 
lows. L + is the intersection of the line E\A with the 
sphere. L_, £3, E2, C2 and C3, respectively, are 
the intersections of E\A! , E\B, E X C \ ~E\D, E 1 E 
and EyO with the sphere. We denote the center 
of the sphere by N. We can easily calculate that 
the angle IE1NE3 is 27r/3. Other angles can be 
calculated. 

Let us give this sphere the coordinates (A, 9), the 
longitude and latitude. The domains of the defini- 
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tion are [0, 2ir) and [0, ±7r/2]. The origin of 9, the 
equator, is the great circle in the xz-plane. The 
north pole is at L + . The origin of A is the great 
circle passing through E\ and L + , and A = at E\ . 

Let us take any point P(x,y,0) in the shape 
plane and connect P and E\ with a straight line. 
Let Pi(xi, yi,z\) be the unique intersection of the 
line and the sphere (6) other than E\. Introducing 
r = x 2 + ?/ 2 , we have 



V 
3r 



Xl 



4r 2 + 3 



cosy, yi 



3r 



4r 2 
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4r 2 + 3 



and 



tan# = — . 

x 

If we represent the {x\, yi, z\) in terms of (A, 
have 

^ ■ x a ^ ■ x ■ a 

X} = — sin A cosy, wi = — sin A sin ft, 

zi = -^-(1 + cos A). 
From Eqs. (J7J) and ©, we get 
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we 



sin A 



4V3V^" 



4(a; 2 



y 



(9) 



(10) 



Thus, Eq. (|8J) and Eq. (|10[1 represent the transfor- 
mation between (x,y) and (X,9). 

In the above, we have moved from the shape 
plane to the shape sphere in the configuration 
space. Now, in the momentum space, let us move 
from the shape plane to the shape sphere. We need 
some preparations. We know that three momen- 
tum vectors make a triangle by Eq. (0). We call 
this triangle the momentum triangle. We ex- 
press this triangle in the (£, 77) plane where the £- 
and 77-axes correspond, respectively, to the x- and 
y-axes in Fig. 1. We borrow the notation from Fig. 
1. 

We normalize the length of P2 so that | p>2 1 = 1- 
Further, we rotate P2 so that it aligns with the 
£-axis. Finally, we put the starting point of vec- 
tor P2 at B(— 0.5, 0) and the end-point at C(0.5, 0) 
(see Fig. 1). We denote this vector by p>2. Corre- 
spondingly, P3 and pi arc transformed to P3 and 
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Pi. Then, P3 starts at the end-point of p>2, Pi 
starts at the end-point of p>3, and the end-point of 
Pi returns to B. 

pi = PS, p 2 = B3, p 3 = UP. (11) 

The (£, 77) plane can be called the shape plane of 
momentum triangle. On this plane, we put the 
sphere of radius V3/4 as in Fig. 1. We intro- 
duce the longitude A and latitude 6 on this sphere. 
Then the transition from the shape plane (£, rf) to 
the shape sphere (A, 0) in the momentum space 
can be carried out perfectly similar to the case of 
the configuration space. We obtain the transforma- 
tion equations jS]) and (fTTJ|) with (A, O) instead of 
(A, 6) and with (£,77) instead of (x,y). 

In the next subsection we give this triangle the 
size and orientation. 

3.4 The remaining variables 

In subsection 3.2, we obtained a set of variables 
(f , F, R*, f2*). In section 3.3, we expressed f = 
(A, 0) and F = (A, 9). In this subsection, we 
change R* to a more convenient variable. Wc 
briefly talk about f2* in the last part of this sec- 
tion. 

R* represents the relative size of configuration 
and momentum triangles. This is obviously related 
to the relative magnitude of the (absolute value of) 
potential energy U and the kinetic energy T (see 
Eqs. (J5J) and ©). For a given energy and a given 
configuration of three bodies, the configuration tri- 
angle is smaller if U is larger, whereas, the mo- 
mentum triangle is smaller if T is smaller. This 
consideration makes it plausible to use virial ratio 
k to parametrize the relative size of configuration 
and momentum triangles, k is defined by 

T 

k=-. (12) 
U 1 

There are two advantages in using the virial ratio 
as one of the variables. One is that the global prop- 
erty of the system can be easily grasped. In fact, 
the total energy h of the system is positive if k > 1, 
whereas h is negative if k < 1. The other advan- 
tage is that any triple system with negative (resp. 
positive) total energy can be brought to a system 
with a fixed negative (resp. positive) energy by a 
similarity transformation. Indeed, we used in §3.2 



a scale transformation when we normalize the size 
of triangles. In that transformation, it turns out 
/3 = or 1 ! 2 . This transformation is equivalent with 
that of the total energy. Let us use this fact to 
connect k and R*. 

Following the notation of §3.2 and §3.3, we have 
r = |q 2 — qa| and R = |p 2 |. When we change the 
scale r,R^r 1,R*, then the relation between these 
variables is R* = y/rR. The relation between k and 
r, R, R* is 

T(fl,A,9) _ T(R*,A,Q) 
U(r,X,9) U(1,X,6) ' 1 6) 

where T and U are represented as functions of 
(R,A,Q) and (r,A,0). 

We not yet determine the relative angle of the 
x-axis and £-axis. This is related to the starting 
position of measuring angle Q* between configura- 
tion and momentum triangles. We take the £-axis 
(i.e., P2) along q 2 . ft* is then the angle between 
P2 and P2 or the angle between q 2 and p 2 , so 

wr K ^wr Kin '^\ (14) 

where lZ(fl) is the rotational matrix around the z- 
axis by fi. 

Finally, we have six variables for the planer three- 
body problem (A, 6, A, ©, k, oS) where we use u> in- 
stead of n*. This space is S 2 x S 2 x (I x S 1 ) where 
I = [0, 1) or / = [0, 00]. If we consider systems with 
negative energy, then I = [0, 1) and all variables are 
bounded. 

3.5 Example orbits in the shape 
space 

Let us represent some preiodic orbits in our shape 
space. 

Collincar motions are on the equators of the con- 
figuration and momentum shape spheres, and u is 
equal to or 7r, whereas k changes with time. We 
know that the isosceles motions move on the merid- 
ians of the configuration shape sphere. These mo- 
tions are represented also on the meridians of the 
momentum shape sphere. 

The position (A, 9) of the Euler collinear motion 
is fixed on the equator of the configuration shape 
and changes depending on the mass ratio of bodies. 
Its position (A, O) is also fixed depending on the 
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mass ratio and (A, 9), and generally 0^0. The 
Lagrange motion is fixed at the north or south pole 
of the configuration shape sphere (See th ecross in 
Fig. 2). (A, 0) is fixed and depends only by the 
mass ratio. Generally 0/0. See Fig. 3 for the 
motion of (k, oS). The loci are concentric closed 
curves with center at k = 0.5, u = tt/2 which is the 
position of the so-called Lagrange solution. 

In the case of the famous figure-eight solution, 
the motion is represnted by the identical curves 
in the X-9 sphere and A-0 sphere because of the 
similarity of configuration and momentum triangles 
j'20j . We show this orbit on the ccy-plane in Fig. 
Its k-u motion is interesting. We show these in 
the Fig. [3l Blue and red triangles correspond to 
collincar configuration, whereas a green triangle to 
isosceles configuration. 

In the above three examples, it is interesting to 
note that these orbits respectively have the same 
trajectories in the configuration and momentum 
spheres. The difference manifest itself in the (k, w)- 
surface. In a sense, these orbits are degenerate. We 
expect that general periodic orbits have different 
trajectories in three surfaces of the shape space. 
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Figure 2: The figure for the orbit of figure-eight on 
the xy-planc. 
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Figure 3: The figure for the orbit of figure-eight 
and lagrange with k=0.1, 0.2, 0.5 on the kw-planc. 

4 The final step to our vari- 
ables 

The space represented by variables (A, 9, A, 0, k, u>) 
can be interpreted in two different ways. The first 
interpretation is that this is a shape space. This is 
an extension of the shape sphere ([T3], [7]). The 
other interpretation is that this is an initial con- 
dition space. Any possible initial condition of the 
planar three-body problem can be expressed in this 
space. 



In this section, we try to further decrease the num- 
ber of variables by one. This can be realized if we 
find a global surface of section. To find a global 
surface of section is motivated from the rectilin- 
ear problem ([ID])- We suggest the surface 1 = 0, 
where / is the moment of inertia of the triple sys- 
tem and / is its time derivative. 

We may change variables from (A, 9, A, 0, k, w) 
to (A, 9, L, I, k, lo) where L is the total angular mo- 
mentum of the system. The equations for the trans- 
formation are 

l = <n A (ftMpO = R * J2 * A c^MPi) ( 15 ) 

i i 

/ = ^> ■ (K(u)pi) = R*J2*i- (ftMPi), (16) 

i i 



4.1 Global surface of section 



G 



where 1Z(lu) is the rotation matrix used in §3.4. 

Now we state the following result. Then we 
can move to the hyper-surface (A, 0, L, I = 0, k, ui). 
This hyper-surface plays a role of the global surface 
of section in the sense that almost all (in the sense 
of measure) orbits pass through this surface. 

Proposition. Orbits except those of measure zero 
experience 1 = 0. 

Proof. Let us first recall the classification of [6]: 
H; hyperbolic motions of three bodies, 
HE^; hyperbolic-elliptic motions 
in which particle i escape, 
P; parabolic motions of three bodies, 
PEj-; parabolic-elliptic motions 
in which particle j escape, 
B; bounded motions, 
OS; oscillatory motions. 

These motions can be initial motions and final 
motions. Then there are 36 combinations of ini- 
tial and final motions such as H~ - B + where '— ' 
indicates the initial motion and '+' indicates the 
final motion. Among these, combinations of escape 
initial motions and escape final motions give nec- 
essarily 1 = because I — ► oo as £ — > ±oo and I 
should attain a minimum at some finite t. In the 
case of OS, there is a minimum of I in each oscil- 
lation, so there is a time such that 1 = 0. 

The remaining combinations to be examined are 
between B and escape motions and between B it- 
self. However, the combinations of B and escape 
motions occupy a set of measure zero in the phase 
space, as has been shown ([2]). Therefore we only 
need to check the last combination, B~ - B + . 

If the orbit does not experience I = in the 
range to < t < oo, then it means / > or / < 
in this range of t. We put J sup = sup to<t<00 / and 
7j n f = inf to<t<00 I. Then for the case I > 0, 

I(t) = f i(t)dt > I(t )+i inl (t-t Q ), (17) 

and for the case I < 0, 

/(f) = f Ht)dt< I(t ) -|/ 8up |(f-fo)(18) 

Let us first consider the case (fTT)) . In order that 
I(t) remains finite as t — - > oo, it is necessary that 



hni = 0. This means that the triple system asymp- 
totically tends to a configuration with constant I 
as t — > oo. In order that I(t) remains greater than 
zero in the case (|T8"| as t — > oo, it is necessary that 
/sup = 0. In this case, the triple system asymptot- 
ically tends to a configuration with constant / as 
t — > oo. If I(t) — > as t increases, this corresponds 
to triple collision with a finite collision time t* . In 
this case, the system does not experience I = for 
to < t < t* . As is well known, orbits which expe- 
rience triple collision occupy a set of zero-measure 
in the phase space. 

The remaining problem is to estimate the mea- 
sure of the orbits which asymptotically tend to 
I = as t — * oo. These orbits constitute the stable 
set of orbits on the set 1 = 0. In any situation, this 
stable set does not span the phase space since oth- 
erwise contradiction to the volume preservation will 
be derived. Therefore, the set of orbits which does 
not experience I = is of measure zero. Q.E.D. 

4.2 An extension to the three- 
dimensional case 

Before extending our results to the three dimen- 
sions, we discuss on the dimension of the initial 
value space of the three-dimensional three-body 
problem. There are ten equivalence relations. Six 
of these are for the translations, three of these are 
for the rotation, and one of these are for the scale. 
Each of them has its equivalence class. We denote 
these sets by E^(i = 1,2,. ..,10), dim Ej = 1 and 
Ej n Ej = 0. We define the set X as R 18 /(LI I 1 £ 1 E i ), 
then the dimension of X is 8. We will find eight 
bounded variables in X in what follows. 

Even in the three-dimensional case, we can de- 
fine the plane in which the three bodies live. We 
take this plane as the initial configuration plane. 
However, in contrast to the planar case, the plane 
in which the configuration triangle exists and the 
plane in which the momentum triangle exists moves 
separately. Thus we need to specify the relative po- 
sition of two planes. 

Now, we look for the remaining variables. Let us 
consider an arbitrary state of the triple system, and 
suppose that the shapes and the sizes of configura- 
tion and momentum triangles are determined. As 
before, we put the center of gravity of the configu- 
ration triangle at the origin of the (x, y)-plane, and 
put q3 — q2 along the x-axis. We here define the 
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nominal position of the momentum triangle. Let 
Pi be the momentum of particle i in the nominal 
position. We put the center of mass of the triangle 
made with at the origin of the (x, y)-plane, and 
put P2 parallel to q2. We take the £-axis in this 
direction, and take the ry-axis perpendicular to the 
£-axis. The £-axis is defined as the third coordinate 
axis in the momentum space. We denote the (£, 77)- 
planc by ttq. So, initally, the plane of the nominal 
momentum triangle coincides with the plane of the 
configuration triangle. 

In order to bring the nominal momentum triangle 
to the the actual momentum triangle, we need three 
steps. 

(1) Rotate 7To around the £-axis by angle w. Wc 
call the new plane tt' , and denote the £-axis of 
tt' by £v 

(2) Rotate tt' around the £ n i by angle ip. We call 
the new plane tt" . Denote the £-axis of tt" by 

(3) Rotate tt" around C w » by angle </>. Then the 
triangle coincides with the actual momentum 
triangle. 

Here, the angle u> coincides with that in the two- 
dimensional case. 

Now it is easy to make the initial conditions 
for the three-dimensional case. At first, we 
make the two-dimensional initial condition for the 
form of the momentum triangle at the nominal 
place. Second, we rotate the plane to the posi- 
tion where the actual momentum triangle exists. 
Then the eight variables for the three-dimensional 
case arc (A,9,A,Q,k,u;,ip,<p), As in the two di- 
mensional case we can transform variables from 
(A, 9, A, 0, k, ui, ip, </>) to (A, 9, L, I, k, u>, ip, <p), where 
L and / are the total angular momentum and the 
derivative of the total moment of inertia. Finally, 
the surface defined by the 1 = will be the global 
surface of section as in the two-dimensional case. 

5 Summary and Discussion 

In this report, for a given mass, we have extended 
the free-fall problem (FFP) to the full three-body 
problem of three dimensions. We find new vari- 
ables (A, 9, L, I, k, u>, (j>, ip) for the three-dimensional 



case and (X,9,L,I,k,uj) for the two-dimensional 
case which are convenient for computer simula- 
tions. The reasons for it are the following. 

(1) If the virial ratio k is positive and large, then 
the total energy of the triple system is positive 
and its final motion is simple. So we can omit 
this case from our consideration, and then we 
can bound the value of k. 

(2) If we set the virial ratio k < 1, then the domain 
of the definition for L and I is bounded. 

(3) The surface defined by / = in the space 
(A, 9, L, J, k, u), <f>, ip) becomes the global sur- 
face of section. We can map the structure of 
the whole phase space on this surface of sec- 
tion. 

Let us briefly discuss about the future role of the 
present result. The progress of computer enables 
us to calculate immense number of orbits. We can 
see the projections of the phase space through in- 
tegrations of orbits by fixing some of the variables 
such as k, L, and / in the initial value space. 

As seen in §3.5, known periodic orbits occupy 
special positions in the shape space. We hope that 
new kinds of periodic orbits may be found with the 
aid of this shape-space representation. 
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